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Abstract 

We study high-momentum distributions and short-range correlation probabilities in the deuteron 
with a variety of modern potentials based on chiral effective field theory up to fifth order in the 
chiral expansion. We also consider some conventional (meson-exchange and phenomenological) 
interactions. We examine our predictions in the context of short-range correlation probabilities as 
extracted from analyses of inclusive electron scattering data and discuss whether modern interac¬ 
tions can be reconciled with the latter. 
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I. INTRODUCTION 


High-momentum components in the nuclear wave function and in momentum distribu¬ 
tions are a reflection of short-range correlations (SRC) in nuclei. The presence of high- 
momentum components is mainly due to the repulsive short-range central force and the 
tensor force. 

Although high-quality meson-theoretic interactions continue to be employed in contem¬ 
porary calculations of nuclear structure and reactions, since the 1990’s there exists a gen¬ 
eral understanding that chiral effective held theory (EFT) PQ [2] is a superior framework. 
First, chiral EFT has a hrm connection with quantum chromodynamics (QCD) through the 
symmetries of low-energy QCD. Second, it allows for a systematic expansion which makes 
possible a quantihcation of the theoretical error. At each order of chiral perturbation theory 
(yPT), the uncertainty associated with a particular prediction can be controlled and quan- 
tihed. For these reasons, nuclear chiral effective theory is becoming increasingly popular as 
a model-independent approach. 

The NN potentials constructed within chiral EFT are generally softer than “conven¬ 
tional” potentials, which makes them computationally more amenable to nuclear structure 
calculations. Also, potentials with a low resolution scale, obtained through a unitary trans¬ 
formation (that is, through renormalization group (RG) methods [3]) applied to a “harder” 
interaction, are very popular for many-body calculations. The resulting “low-k” potentials 
are equivalent to the original ones for all physical purposes, although essentially void of high 
momentum components. 

On the experimental side, inclusive electron scattering measurements at high momentum 
transfer, on both light and heavy nuclei, have been analyzed with the purpose of extracting 
information on short-range correlations Hi]. In a suitable range of and xb , the cross 
section is factorized in order to single out the probability of a nucleon to be involved in 
SRC, either two-body or three-body. When extended to nuclear matter, this probability 
is equivalent to the “wound integral”, which measures the amount of correlations in the 
so-called defect function [7|. Information about two-body correlations can also be obtained 
in coincidence experiments involving knock-out of a nucleon pair with protons [8] or elec¬ 
trons [9HI2|. 

Nuclear scaling and the plateaus seen in inclusive scattering cross section ratios H are 
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due to the dominance of SRC for momenta above approximately 2 fm“^. In the same region, 
the momentum distribution in a nucleus relative to the one in the deuteron becomes almost 
flat, so that those distributions simply scale with A. 

The discussion around some of these measurements which is presently going on in the 
literature is quite intriguing. The probabilities mentioned above are a manifestation of 
the off-shell nature of the potential, which cannot be determined uniquely from NN elastic 
data and is not an observable. Interactions may differ dramatically in their off-shell behavior 
while remaining phase-equivalent. The most striking example is provided by the RG-evolved 
potentials we mentioned above, where the high-momentum structures of the original and 
the RG-evolved potential are obviously not the same. 

Naturally, a low-resolution scale will impact the ability to resolve high momentum regions. 
With regard to this point, it has been noted that, if a unitary transformation is applied to 
both wave functions and operators, one regains the invariance of the cross section, as one 
should [ 13 ], thus attaining a consistent description of short- and long-range physics. This 
has been addressed recently by Neff et al., who show how the short-range information can 
be recovered by transforming the density operators m. 

On the other hand, chiral potentials such as those developed in Refs. USHIZI are not 
low-momentum interactions in the sense of a Viow-k- In this paper, we examine those from 
the point of view of SRC. To broaden the discussion, we start with an analysis of SRC and 
conventional (that is, non-chiral) NN potentials, including high-precision potentials from the 
1990’s as well as a phenomenological one. We then move to a similar analysis with chiral 
interactions. 

We take the deuteron as our sample system. We recall that the high-momentum part 
of the momentum distribution shows similar features in nuclei with ^4=2 to 40 [18]. Thus, 
the deuteron offers representative features. Furthermore, deuteron SRC probabilities are 
a crucial element in the estimation of SRC probabilities in heavier nuclei as obtained in 

Ref. [ 3 ]. 

We calculate the momentum distribution and the probability of SRC up to 5th order of 
chiral effective theory. Working with the ^4=2 system, we can go to any order of chiral EFT 
where NN potentials are available, without the need to worry about the corresponding three- 
nucleon forces (3NF). Although 3NF at N^LO and N^LO have been worked out [TUIET] . their 
application in few- and many-body systems still presents considerable challenges and requires 
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TABLE I: Probabilities of SRC and D-state probabilities for the potentials considered in Fig. 1. 


Model 

a2N{d) 

Pd 

CD-Bonn 

0.032 

0.0485 

AV18 

0.042 

0.0578 

Nijmegen II 

0.041 

0.0565 


unavoidable omissions/approximations. Calculations in the deuteron are free of those and 
thus well-controlled. We address cutoff dependence and order-by-order convergence from 
lowest to 5th order of the chiral expansion. 

Some of the questions we wish to address are: To which extent are modern, non- 
phenomenological interactions (chiral or not) consistent with the information as extracted 
from A{e, e')X measurements? What do we learn, on fundamental grounds, from the answer 
to this question? Are there characteristic differences among particular families of potentials 
from which we can obtain physical insight (beyond phenomenological observations)? 

Our hndings and conclusions are summarized in Section 


II. HIGH-MOMENTUM DISTRIBUTION IN THE DEUTERON 


A. Meson theory and phenomenology 


We begin with a step back into the 1990’s by considering three members of the “high- 
precision” family of NN potentials, namely CD-Bonn [ 22 ], Nijmegen II [ 22 ], and AV18 [ 21 ]. 
The respective momentum distributions in the deuteron, with focus on high-momentum 
components, are shown in Fig. 1. p{k) is the Fourier transform squared of the coordinate 
space wave function. There are noticeable differences between the (softer) predictions from 
CD-Bonn and those from the other two potentials, which are essentially indistinguishable. 

To see how these differences carry into the probability of SRC, we follow Ref. [18] and 
dehne the probability of SRC in the deuteron as 


a2N{d) = in 



p{k)k‘^ dk , 


( 1 ) 
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FIG. 1: (Color online) High momentum distributions in the deuteron as predicted by: CD-Bonn 
(solid red); Nijmegen II (dotted blue); and AV18 (dashed green). 

where kmin is taken to be 1.4 fm“^ (276 MeV). This definition was adopted in Ref. |1], 
where the choice of the lower integration limit is suggested by the onset of scaling of the 
cross section, which should signal the dominance of scattering from a strongly correlated 
nucleon. In Ref. [ 1 ], the ratio of the per-nucleon probability of two nucleon (2N) SRC in a 
nucleus relative to ^He is argued to be equal to the ratio of the inclusive electron scattering 
cross sections in the appropriate scaling region. The absolute per-nucleon probability in 
a nucleus can then be deduced if the absolute per-nucleon probability in ^He is known. 
The latter is the product of the absolute per-nucleon probability in the deuteron, stated 
as 0.041T 0.008 in Ref. [25], and the relative probability of 2N SRC in ^He relative to the 
deuteron. Namely, 

fl 2 Af(^) = cb 2 N{-^/^He)a 2 N^He) and a 2 N^He) = a 2 N^He/(£)a 2 N{d) ■ (2) 

We note that the values for the deuteron and the ratio of ^He to deuteron contain some 
model dependence from theoretical calculations (see Ref. |1| and Refs. [2,6,15,16] therein), 
likely to propagate in the predictions for heavier nuclei. 


5 



In Table 1^ we show the probability as dehned in Eq. Q for the interactions used in Fig. 1. 
As an additional, related information, we also show the corresponding D-state probabilities. 
As expected in light of Fig. there is a signihcant difference between CD-Bonn and the 
other two cases, with the AV18 and Nijmegen II predictions closer to the value used in the 
analysis from Ref. [25]. 

The differences noted above are due to the non-local nature of CD-Bonn, which adopts 
fully relativistic momentum-space expressions for the one-pion-exchange. More precisely, the 
off-shell nature of CD-Bonn is based upon the relativistic Feynman amplitudes for meson 
exchange. This determines well-founded non-localities in the tensor force, whereas Nijmegen 
II and AV18 make use of the non-relativistic, static one-pion-exchange which generates a 
local tensor force. The characteristically softer nature of a relativistic momentum-space 
potential reflected in Table [I] is a desirable feature for the purpose of applications in nuclear 
structure. 

We will come back to this point for a more complete discussion after the next section. 

B. Interactions based on chiral EFT 

In spite of the good theoretical foundation behind meson-exchange Feynman amplitudes, 
meson theory does not provide a systematic approach to constructing nuclear forces. As 
mentioned in the Introduction, chiral EFT presents the opportunity for such systematic 
development. 

Crucial for a nuclear EFT are the processes of regularization and renormalization. Con¬ 
cerning the former, all chiral interactions are multiplied by a regulator function which typi¬ 
cally has the form; 

/(p',p)=exp[-(p7A)2--(p/A)2T, (3) 

where A is known as the cutoff parameter. 

Nucleon-nucleon potentials have been developed at different orders and cutoff values [161 
HU. Chiral EFT predictions allow for the quantihcation of uncertainties that stem from the 
truncation error and cutoff variations (as well as additional sources of errors). 

Consistent with that philosophy, in Fig. we show the momentum distribution in the 
deuteron, including Eve orders of the chiral expansion. We note that the potential at N^LO 
is a preliminary version of a high-precision nucleon-nucleon potential at fifth order EH- 
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FIG. 2: (Color online) Momentum distributions in the deuteron predicted with chiral potentials 
at: LO (dotted); NLO (dash-double dot); N^LO (dash-dot); N^LO (dash); and N^LO (solid). The 
cutoff is fixed at 500 MeV. 

(Some wavy structures noticeable in the figure are most likely due to the polynomial nature 
of the EFT contacts.) We observe huge variations at the lowest orders, particularly from 
LO to NLO, and a clear convergence pattern with increasing order. 

Table |IT| shows the integrated probabilities corresponding to Fig. At the higher, con¬ 
verging orders, the SRC probabilities are not very different from the one predicted by CD- 
Bonn. That is, chiral potentials with suitable cutoff can be constructed with excellent fit 
to the NN data and off-shell nature similar to the highest-quality non-local meson-exchange 
forces. 

The truncation error at order n is defined as the difference between the predictions at 
orders n + 1 and n. Thus, the error of a 2 N{d) at N^LO is ±0.004. We have also considered 
variations of the cutoff parameter between 500 and 600 MeV. At N^LO, we found the value 
of a 27 v(d) with cutoff of 600 MeV to be nearly the same as with 500 MeV. Therefore, cutoff 
uncertainty is below the truncation error, and our final result at N^LO can be stated as 
0.030±0.004. Concerning the uncertainty of the N^LO result, the prediction for the next 
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TABLE II: Probabilities of SRC and D-state probabilities for the chiral interactions considered in 

Fig. [3 


Model 

a2N{d) 

Pd 

LO 

0.047 

0.0757 

NLO 

0.015 

0.0313 

N^LO 

0.022 

0.0417 

N^LO 

0.030 

0.0451 

N^LO 

0.026 

0.0414 


higher order is unknown. Therefore, assuming (pessimistically) the same truncation error 
as at N^LO, the prediction at N^LO can be stated as 0.026±0.004. 

C. Discussion 


The deuteron is the simplest system where off-shell behavior can be explored. Charac¬ 
teristic differences exist between meson-theoretic potentials using fully relativistic one-pion 
exchange amplitudes (that is, non-local tensor forces) and those which use static one-pion- 
exchange. Off-shell behavior is not observable and thus cannot be uniquely determined by 
measurements. The best one can do is to have a good theoretical foundation for it. In meson 
theory, this is provided by relativistic meson-exchange amplitudes. 


In Section IIA we saw that the SRC probability predicted with CD-Bonn is roughly 25% 
below the value of a 2 N{d) cited in Ref. |1] and used to evaluate absolute probabilities in 
heavier nuclei. A discrepancy of qualitatively similar nature exists for the wound integral 
in nuclear matter. More precisely, conventional non-local potentials are known to predict 
about 10-15% for the wound integral (at normal density) [7], whereas a value of about 25% 
is cited from extrapolation to nuclear matter of the empirical information [28]. In fact, the 
ratio of the probability for nucleus A to the deuteron, extrapolated to inhnite symmetric 








matter, is given to be 6.5±1.0 [28]. With the absolute probability for the deuteron taken 
equal to 0.04, the value cited above is obtained. However, when the a2Ar((i)=0.032 value 
from CD-Bonn is used, one obtains 20%, and with a2Ar(d)=0.026 as from the converged 
chiral results, a value of 17% is obtained for empirical short-range correlations in nuclear 
matter, which is getting closer to the predicted wound integrals. 

In summary, the question of consistency between description of short and long range 
physics seems to go beyond the (intrinsically) “low-momentum” nature of some potentials. 
Instead, it points to non-locality in the tensor force, a feature which has been found since a 
long time to be very attractive in nuclear structure. So, we are revisiting an old issue which 
resurfaces in the light of new experiments. We suggest to take a broad view of it, combining 
our former and present understanding of microscopic nuclear forces and their development. 


III. CONCLUSIONS AND OUTLOOK 

The deuteron is a beautifully simple benchmark for theories of nuclear forces. From this 
study, we conclude that predictions of high-momentum distributions in the deuteron with 
high-quality non-local meson-exchange forces or state-of-the-art chiral forces are systemat¬ 
ically lower than what is used to extract empirical information for heavier nuclei. Taking 
those results into account leads to a better agreement between SRC in nuclear matter and 
theoretical predictions. 

We plan to extend our microscopic analysis to the A=3 system using a broad spectrum 
of interactions as in the present study. We hope this will shed more light on how to reconcile 
theory and empirical analyses. 

Finally, some caution should be excercised in the interpretation of the empirical informa¬ 
tion discussed above as an experimental constraint on the off-shell behavior. 
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